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ABSTRACT 

The standard hght curve of a microlensing event provides only two constraints 
on the six unknown parameters of the lens. We show that narrow-band 
photometry during a microlensing event of a giant star can in addition determine 
the angular radius of the Einstein ring and the proper motion of the lens. This 
possibility results from the fact that the extended atmospheres of giants emit 
the cores of resonant lines primarily from a narrow ring (limb brightening). 
The radius of the emission ring can be determined to a precision of ±20% by 
detailed spectroscopic observations of the source after the lensing event has 
ended. A considerable fraction of the clump giants in the bulge have a ring 
radius ^ 10^^ cm, within the range of Einstein radii for sub-solar mass lenses. 
The extended thin ring also provides a sensitive probe of possible planetary 
companions as it sweeps across the lens plane. The ring signature can be detected 
photometrically, using a narrow-band filter centered on the Call K line at 3933A. 
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1. INTRODUCTION 



By now there are more than 50 microlensing events detected towards the Large 
Magellanic Cloud (LMC) and the Galactic bulge (Alcock et al. 1995; Auborg et al. 1993; 
Udalski et al. 1994a). The standard light curve of an event is characterized by two 
observables: its peak amplitude and duration (Paczyhski 1986). However, there are six 
unknown parameters that define the physical properties of the lens, namely its mass, distance, 
its three velocity components, and the impact parameter of its trajectory relative to the line 
of sight to the source. The mass, distance, and transverse velocity of the lens combine 
to set the duration of the event, and together with the impact parameter they also define 
its peak amplitude. Breaking the degeneracy among these parameters is of fundamental 
importance to the interpretation of microlensing light curves. For observations towards the 
LMC, it is essential to discriminate between lenses that are located within the LMC (Sahu 
1994) and lenses which are part of the Galactic halo or thick disk (Gould, Miralda-Escude & 
Bahcall 1994); for bulge observations, it is important to separate lenses that are embedded 
in the Galactic disk from lenses inside a possible Galactic bar (Zhao, Spergel, & Rich 1995; 
Han & Gould 1995). Removal of the existing degeneracy between the lens parameters is 
a necessary step towards a better understanding of the nature of the lenses. A variety of 
methods to partially achieve this goal were suggested recently. The most ambitious proposal 
involves a parallax satellite mission to measure the lens proper motion (Refsdal 1966; Gould 
1992, 1994b, 1995a). Other proposals consider alterations of the standard light curve due 
to different properties of the source, such as its finite size (Gould 1994a, 1995b; Nemiroff 
& Wickramasinghe 1994; Witt & Mao 1994), its polarized emission (Simmons, Newsam, 
& Willis 1994), spectral shifts due to its rotation and finite size (Maoz & Gould 1994), or 
shifts in its color due to a luminous lens (Kamionkowski 1995; Udalski et al. 1994c). In this 
Letter, we show that narrow-band photometry of lensed giant stars can be used to measure 
the proper motion and the angular radius of the Einstein ring of the lens. 

Optical starlight is dominated by photospheric continuum emission. The continuum 
emission generally peaks at the projected center of the stellar disk and falls-off towards the 
stellar limb. In cool stars, some spectral lines (such as Ha or Ca II) are still optically thick 
at the limb and therefore appear strongly in emission there. This resonant line scattering 
was predicted for the Sun in 1870 by C.A. Young and was later observed during solar eclipses 
(Zirin 1988). The limb brightening of various strong lines in the Sun extends only over a 
ring of width ^ O.Oli?©. Cool stars with a much smaller surface gravity, such as giants 
and supergiants, have atmospheres which are more extended than the Sun by 1-2 orders of 
magnitude and therefore show a much stronger resonant line scattering. Limb brightening 
of supergiants has been observed directly only recently (Hannif et al. 1992). 
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For giant stars, the line emission ring has a radius of 0.05-lAU, comparable to the 
Einstein radius of sub-solar mass O.IMq) lenses. A ring of this size provides just the 
baseline necessary to mimic observations of a microlensing parallax by an Earth-satellite 
system. Eor a specific star, the actual radius of the ring can be measured up to a precision 
of ~ 20% after the lensing event has ended, using traditional methods of high-resolution 
spectroscopy and infrared photometry (see McWilliam & Rich 1994 for application to 
Galactic bulge giants). Based on the measured radius, it is then possible to determine 
the projected Einstein radius and the proper motion of the lens relative to the line of sight. 
The extended thin emission ring also provides a sensitive probe of planetary companions as 
it sweeps across the lens plane. 

The outline of this paper is as follows. In §2 we present the method used for calculating 
the narrow-band light curves of giants. We calculate disk intensity distributions for several 
atmospheric models and a 6-level Ca II atom model. In §3 we estimate the probability 
distribution of radii for clump giant stars in the Galactic bulge, based on the color-magnitude 
diagrams of the OGLE observations (Stanek et al. 1994). Numerical examples for narrow- 
band light curves of bulge giants are presented in §4. Einally, §5 summarizes the main 
conclusions of this work. 



2. NARROW-BAND LIGHT CURVES OF GIANTS 



The amplification of a point source by a point mass, Af , depends only on their projected 
separation d (Paczyhski 1986), 

where d is expressed in units of the angular radius of the Einstein ring of the lens, 
9^ = [[4:GM/c^][D^g/ D^^D^g])^^'^ ^ and D^^^D^^^ and D^^ are the distances between the 
observer, lens and source. The total fiux received from an extended source is therefore 
obtained by integration over its infinitesimal elements, 

F{t)= r\drB{r) d0A{d:), (2) 
Jo Jo 

where B[r) is the surface brightness profile of the source in the projected polar coordinates 
(r, 6) around its center. Eor simplicity, we assume that the stellar emission is circularly 
symmetric, and denote the source radius by Rs = Rs X (O^D^g). All projected length scales 
are normalized by the Einstein ring radius. Using the law of cosines 

d = \do -f\ = {dl - 2dor cos e + r' fl^ , = + v'efl^, (3) 



- 4 - 



where do is the projected separation between the lens and the source center, b is the projected 
impact parameter of the source center, v is the transverse velocity of the lens in units of 
^_B-^oL5 t = is the midpoint time of the lensing event. For a circularly symmetric 
source, the resulting light curve is time-symmetric with F[t) = F[—t). 




Fig. 1. — Stellar disk intensity distributions from a model atmosphere calculated in spherical 
geometry. The solid line shows the center-to-limb intensity variation in the core of the Ca II K 
A 3933 line for a 100 Rq cool (4150 K) giant star. The dotted line shows the broad-band 
continuum intensity variation in the spectral region adjacent to the Ca II line. 

The narrow-band light curve of a giant is obtained by substituting its center-to-limb 
intensity distribution, i?(r), in equation (2). In calculating B[r) we ignore any surface 
features on the otherwise circularly-symmetric stellar disk, i.e. our radiative transfer 
models are one-dimensional. As initial conditions for constructing spherical atmospheres 
we use constant-flux plane-parallel blanketed photospheres from Kurucz (1993), and outer 
atmosphere models from Dupree, Hartmann, & Avrett (1984) and Harper (1992). The multi 
level non-LTF radiative transfer in these models is then solved in a spherical geometry using 
the SMULTI code (Harper 1994). We performed calculations with several atomic models 
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which have strong observable resonance hnes (such as H, Ca II, Mg II) and found that both 
the Ca II H & K hnes show significant scattering above the hmb, and are idealy suited for 
photometry. The model used for the Ca II atom has 6 levels; it treats both resonance lines 
and the infrared triplet in detail. A discussion of a similar calculation and a comparison 
with plane-parallel profiles can be found in Sasselov & Karovska (1994). An example for the 
most extended atmosphere we have considered is illustrated in Fig. 1, where we show the 
intensity distribution across a giant with a radius of lOOi?© and an effective temperature of 
4150 K. The intensity distribution for the Ca II K line core refiects our optimal choice of a 
wavelength band. This filter passband was chosen based on our modeled line profiles and on 
practical considerations regarding existing observational setups-to be discussed in § 4. 

3. GIANTS IN THE BULGE 

The giants in the bulge have a relatively high metalicity (McWilliam & Rich 1994), 
implying a high opacity in their interiors. For a given stellar mass, a higher opacity favors 
convection and leads to a larger radius, i?^, and a cooler effective temperature (Iben & 
Faulkner 1968). The red giants in the bulge are similar to the members of the old open 
cluster M 67, whose color-magnitude relation is often used to represent the bulge (Ibata 
& Gilmore 1995). In order to determine the probability distribution of giant radii in the 
bulge, we use the color-magnitude diagrams provided by the OGLF collaboration for Baade's 
Window (Stanek et al. 1994; see similar analysis by Witt 1995). 

The relation between absolute magnitude, color, and radius, can be established using 
direct measurements of giant radii in the solar neighborhood and through the calibration of 
the M 67 cluster. In particular, one needs to find an adequate temperature calibration for the 
Cousins V-I color used in the OGLF observations of the bulge. We determine the extinction 
based on the better known value for F(B-V) = 0.48 (e.g. McWilliam & Rich 1994), and 
the transformations of Taylor (1986) for K2 III stars. We find F(V-I) = 0.65. This result 
is in good agreement with the HST photometry of Baade's Window by Holtzman et al. 
(1993)-using the new abundance estimates by McWilliam & Rich (1994), as well as with the 
recent result from W UMa systems by Rucinski (1995). Next, we use as most reliable the 
direct temperature calibration of V-K photometry from lunar-occultation angular-diameters 
of field giants by Ridgway et al. (1980) and from Michelson interferometry by Di Benedetto & 
Rabbia (1987) and Mozurkewich et al. (1991). This calibration is transformed to match the 
bulge giants by using the giants in M 67 as surrogates with data from Taylor & Joner (1988) 
and Joner & Taylor (1988). As an independent check, we compared the recent spectroscopic 
temperatures of 11 bulge giants by McWilliam & Rich (1994) to the above calibration of the 
giant's effective temperature and found a very good match, except for three unusual cases. 
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Fig. 2. — Probability distribution of stellar radii in Baade's Window, derived from the 
OGLE data set (Stanek et al. 1994). The clump giants are clearly separated from the Main 
Sequence. The Einstein radius for a lens mass of 10""^ Mq {32Rq) and a lens mass of 0.1 Mq 
(3I6-R0) are marked, assuming D^^ = 8 kpc and D^-^ = 6 kpc. The Main Sequence is only 
qualitatively illustrated-our radius calibration applies to giants and supergiants. 

Using the above temperature calibration and adopting a distance to the Galactic center 
of = 8 kpc, we have computed the radii of all stars for which the above relations are valid 
in the OGLE data set (~ 50% of the entire stellar population). Eigure 2 shows the probability 
distribution of stellar radii for ~ 90% of the OGLE stellar population (V-I > 0.67). Since 
we are interested primarily in cool giants and supergiants, we did not use a separate radius 
calibration for dwarfs; the Main Sequence is shown in Eig. 2 only to illustrate its relative 
fraction in the observed stellar population. We find that giants with radii ^ 20Rq comprise 
13% of the clump giants in Baade's Window. About 19% of all the stars have radii ^ 5Rq. 
The effect of a possible stellar bar on this distribution was found to be small and comparable 
to other uncertainties. However, since a significant fraction of the lenses may reside in such 
a bar (Zhao, Spergel, & Rich 1995), we adopt a distance to the lens of D^-^ = 6 kpc and a 
transverse lens velocity of 100 km s~^ in the numerical examples that follow. 
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4. RESULTS 



The light curves of ring-hke sources separate into two types: those with a single smooth 
peak and those with two sharp spikes. The first type is obtained when the lens is outside 
the ring trajectory, i.e. b ^ i?^, and the second type is obtained when the lens crosses the 
ring, i.e. b < i?^. 

For the extreme first-type curves with Rs <^ 6, it is possible to solve equation (2) 
analytically. The amplification can then be Taylor expanded to second-order, 

A{d) ^ A{do) + A'ido) x{d- do) + ^A"{do) x{d- do)\ (4) 
where its derivatives are obtained from equation (1) 

= "J2(J2+ 4)3/2 ' ^'Vo) = ^3(^2^4)^5/2 ' l^) 

and where equation (3) gives. 



^2 



{d - do) ^ -r cosO + —sin^ ; {d - dof ^ r'^ cos^ . (6) 

2do 

The light curve of a point source, Fo{t), is obtained by using A{do) instead of A{d) in equation 
(2). The fractional deviation from this light curve due to a source with a nonzero size is 
obtained by substituting equations (4)-(6) into equation (2), 

m^^-^^^^ = m{r') . (7) 

-To 

where (r^) = [Iq^ rdr B{r) r^]/[f(^^' rdr B{r)], and where the temporal profile of the excess 
signal is given by, 

dUdl + 2){dl + 4)^ ■ 

2 

For emission from a narrow ring, A(t) = G{t)Rs . Thus, common events with b ~ 0.5 

2 

yield maximum values of A ~ i?^ , so that a 3% accuracy in the photometry of the light 
curve which includes ten points of observations can reveal the signature of a source radius 
Rs ^ 0.1. For our lens-source parameters, this threshold corresponds to a source radius 
^ 2 X 10^^ cm X (M/O.lM0)i/2, which is in the expected range of radii for clump giants in 
the bulge (cf. Fig. 2). Based on a comparison to the exact solution, we find that equations 
(7) and (8) provide a good description to the ring light curve as long as Rs < b. 
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However, the character of the hght curve changes for Rs > 6, when the ring is able to 
cross the lens. In this case, the event contains two pronounced peaks which are centered 
around the two crossing times of the ring through the lens. These crossing times are defined 
by the condition d = \do — Rs\ = 0. The structure of the peaks depends on i?^; a larger 
Rs makes the peaks narrower because the fraction of the ring area that is being affected by 
the lens gets smaller. By fitting the peak structure and separation to a stellar atmosphere 
model, it is possible to find both the impact parameter b and the source size Rs in units of 
the Einstein radius. 

I , , , I , , , , I , , , , I , , , I 
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Fig. 3. — Flux amplification profiles in a narrow passband centered on the Ca II K line for 
lensing of the giant shown in Fig. I by a O.IMq mass lens. The light curves are shown for 
different values of the dimensionless impact parameter, b. The radius of the Ca II ring of 
the giant is Rs = 0.35 in Finstein ring units. We plot the net photon fiux on Farth for the 
spectrophotometry filter described in § 4. 

Figure 3 presents the distribution of narrow-band light curves as a function of the 
impact parameter b for a lens mass M = O.IMq positioned at D^-^ = 6 kpc and moving with 
a transverse velocity of 100 km s~^ relative to the line of sight to a clump giant source with 
Rs = IOOi?0 (cf. Fig. I) at D^^ = 8 kpc. As this figure shows, it is much easier to identify 
the signal of the ring at small impact parameters. Figure 4 compares the narrow-band (solid 
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line) and the broad-band (dashed) hght curves of this source for a lens mass of lO~"^Af0 and 
6 = 0. The main difference between the curves is the existence of two peaks in the narrow- 
band curve. The significance of these peaks is emphasized by the similarity between the 
broad-band curve (dashed) and a standard microlensing curve for a point source with some 
different impact parameter and velocity (dotted). Although the small difference between 
these light curves is detectable, the extraction of the value of Rs from this difference could 
be complicated by observational errors or blending with unlensed starlight in crowded fields 
(Di Stefano & Esin f995). The finite source signature can be measured more reliably and 
unambiguously from the two peaks in its narrow-band light curve. For D^-^ = 6 kpc, the 
fraction of giant lensing events with b < Rg is of order {Rg) 4% X (Af/O.f Af©)"^/^ (cf. Fig. 
2). 



1.4 



1.3 



^ 1 P 
o 

Cl, 

1.1 




Fig. 4. — Comparison of light curves for a lens mass of fO~"^Af0. The narrow-band light 
curve (solid line) is compared to the 4000 A broad-band light curve (dashed), for the giant 
shown in Fig. f with Rs = 3.5 and 6 = 0. The dotted line is a standard microlensing curve 
for a point source, chosen to fit closely the broad-band curve; it corresponds to 6 = f .f5 and 
a transverse lens velocity of 50 km s~^. 



The several decades of experience in fiux measurements of Ca II emission from dwarf 
and giant stars should ease the technical aspects of the narrow-band photometry required 
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here. In practice, one could use a modern version of the four-channel photon-counting 
spectrophotometer {HKP2) of the Mt. Wilson HK project (Duncan et al. 1991) with 20A- 
wide reference bands and 2A-wide H and K bands. For example, we find that the net fiux 
in these H and K passbands per unit area at the stellar surface is 4 X 10^ erg cm~^ s~^, for 
a giant with a radius of 20 i?© and an effective temperature of 4400 K (see Rutten (1984) 
for absolute fiux calibration in these passbands). With the adopted distance to the bulge, 
this translates to a detection rate of about 1 photon/sec on a 2-meter telescope. Therefore, 
high quality (±2%) narrow-band photometry can be achieved using the HKP2 instrument 
with an integration time of ~ 45(-Dt/2 m)~^ minutes on a telescope with a diameter Df 
This result applies to the signal without amplification; the peaks of the light curve shown 
in Fig. 4 are easier to observe. In addition, giants with Rg > 20Rq require much shorter 
integration times as can be seen from Fig. 3. Occasionally, a microlensing event may involve 
a chromospherically active giant which happens to have an active phase at the time of the 
event. Such activity could distort the circular symmetry of the emission ring or add bright 
spots on the stellar disk. However, subsequent photometry would be able to determine 
whether the lensed star was likely to have been active during the event. 

5. CONCLUSIONS 

A considerable fraction of the giants in the bulge have a radius Rs ^ 10^^ cm (cf. Fig. 
2). As illustrated in Fig. 1, the narrow-band (resonance line) images of these giants are thin 
rings. When a ring-like image crosses a lens, its time-symmetric light curve shows two sharp 
peaks. These narrow peaks provide an unambiguous signature of the finite source size. The 
broad-band light curve contains a smooth, weaker signature of the source size (cf. Fig. 4), 
that could be inaccurately interpreted due to observational errors or blending with unlensed 
starlight in crowded fields. The radius of lensed giant can be measured spectroscopically to 
within ±20% after the lensing event has ended. The measured radius can then be combined 
with the narrow-band light curve to infer the proper motion and the angular radius of the 
Finstein ring of the lens. These determinations can be used to discriminate between the 
different distance scales possible for the lens population. 

The existence of an extended thin emission ring enhances the probability of detecting 
low-mass companions of the lens, such as planets (for discussions about planetary lensing 
of a point source, see Mao & Paczyhski 1991, Gould & Loeb 1993, and Bolatto & Falco 
1994). First, the finite source size sweeps through a wide strip in the lens-companion plane. 
Thus, the probability of intersecting the lensing region of infiuence of a companion planet is 
increased relative to the case of a point source by the ratio between the source size Rs and the 

1/2 

Finstein radius of the planet. For bulge observations, this ratio is ~ {Rs/2 X lO^^Afp/ cm). 
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where Mpi is the planet mass in earth mass units. Second, the planetary amplification is 
not averaged-out by the entire stellar disk as in the case of broad-band observations, but 
only by the width of the narrow emission ring-which is an order of magnitude smaller (cf. 
Fig. 4). Even in the extreme case of Fig. 1 where Rs = lOOi?©, the width of the emission 
ring is still smaller than the Finstein radius of a Jupiter-mass planet ~ 2 X 10^^ cm. The 
ring provides two chances for the detection of a planetary-crossing signature. Based on the 
timing of these crossings relative to the primary lens crossings, it is possible to infer both the 
angular separation and the orientation of the companion relative to the primary lens. The 
proper motion of the two objects can be compared for any difference associated with their 
relative orbital motion^. Fach planetary crossing signature may be multiply peaked and is 
expected to last ^ day (cf. Gould & Loeb 1992). Once a lensing event of a giant star has 
triggered one of the currently operating early-warning systems (e.g. Udalski et al. 1994b), 
it would be useful to monitor its light curve on an hourly basis with a narrow-band filter. 
This may prove to be an efficient method in searching for planetary systems around stellar 
lenses. 

We thank R. Di Stefano and S. Mao for useful discussions. 
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